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М.В.ЧУВЕЙКО 


АЛГОРИТМ ВОССТАНОВЛЕНИЯ 
ПЕРИОДИЗИРОВАННОГО СИГНАЛА В РОТОРНЫХ СИСТЕМАХ 


Рассматривается алгоритм восстановления периодической составляющей сигнала в 
системах с аддитивной помехой амплитудного характера при отсутствии априорных 
сведений о частоте восстанавливаемого сигнала. 

Ключевые слова: восстановление сигнала, средняя за период мощность, усилен- 
ный закон больших чисел. 


Введение. Большинство механических систем в своем составе имеют 
узлы, совершающие вращательные движения. При стационарном режиме 
функционирования (неизменные угловые частоты) данных узлов такой вид 
движения характерен тем, что координаты состояния таких систем являют- 
ся периодическими функциями времени или пространственных перемеще- 
ний. Факторы, вызывающие периодические вариации координат состояния, 
разнообразны. При обработке на металлорежущих станках - это разнооб- 
рабатываемость, вызывающая погрешности формы поперечного сечения. 
При рассмотрении качества роторных фрикционных систем - это развитие 
неравномерности износа по периоду вращения. К таким же эффектам при- 
водит дисбаланс роторных систем. Эта особенность позволяет сделать вы- 
вод о том, что проявления возможных дефектов данных узлов также будут 
носить характер, близкий к периодическому. Таким образом, задача диа- 
гностики работоспособности системы может быть сведена к восстановле- 
нию периодических составляющих в выходном сигнале и анализу при- 
частности тех или иных механизмов к их формированию. 

Алгоритм восстановления периодической составляющей. Пусть за- 


дана функция 
Е(@)= Л@+:0, (1) 
где & (1 ) - случайный стационарный (в узком смысле) процесс, удовлетво- 
ряющий свойствам: для № выполняется ЛМ (1) = 0, в (1) и 
(№) независимы при & &; Л(Е)= Л@-Т) - кусочно-моно- 
тонная ограниченная периодическая функция с периодом ТГ. 


Наиболее очевидным и простым методом восстановления формы пе- 
риодической составляющей является операция усреднения с периодом та 


Определим операцию усреднения функции ЁЕ(Г) с периодом 1 
следующим образом: 


1; (®)= бп " Еа+м) 
" о И+1 


При усреднении с периодом Тиз усиленного закона больших чи- 
сел [1] очевидно, что 


(2) 


РУ Г 9 =1. 
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В большинстве случаев точная величина периода 7’ неизвестна, но 
задан некоторый диапазон возможных значений. В таком случае необходи- 
мо провести дополнительное исследование, направленное на получение 
периода усреднения. 


Теорема. Среди множества функций 1; (Г) почти наверное 


наибольшими и равными между собой значениями средней за период мощ- 
НОСТИ 
т 


1 р й 
р(т] = — Л, (14 
Г 0 
обладают те и только те функции, для которых выполняется следующее 
равенство: т = ТИ, где и№ 
Доказательство. 
1. Сначала докажем это для случая, когда Л (Г) = $11(®1+У ), где 


2 И 
@ =—— . Представим функцию следующим образом: /(#)= Пп(е”°). Тогда 
т 
формула (2) примет вид 
п ЦОНу+А) И 
. е . Е (1+ 
9; (1) = ш Пи —_ + ый (3) 
Е: Е " о И+1 
р 211 
где 0 = ——. 
Т 


Исходя из свойств & (Г) можно заключить, что для любого Ё почти 
наверное 








п 
А 
п> (=0 И 
Тогда 
п ве +9) | п ей 
Ош а бб зщ ера 
"> о п+1| п» ной+ | 
Рассмотрим отдельно 
п п-1 
а 2 1+ е“ . 1+ ее“ 
И = Пт = Ни = Пи + ее Пи = еб‘ 
7 цоЙ+[ п+1 "> И+| п п+1 


Но такое возможно лишь при (/= 0, если @ 2МИ. Если же 


@ = 27 И, легко видеть, что 
п 21 
И=Ш —=1. 
"> .0И+ 1 
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Исходя из этого интеграл: 


Ао ь 0.5, 
Е = 
р) т 7.04 | 


если т = Ти; 


если То ТИ. 


2. Теперь докажем это для случая, когда У (Г) - произвольная пе- 


риодическая функция, удовлетворяющая указанным ранее свойствам. 
Представим функцию в виде ряда Фурье: 


Г0= А+ Аз@+у,), 
11 
тогда формула (2) примет вид: 


п 


ОА 


Азш( чу ,). ет 


30 п+1 
10=А+ Ай " эш(фтчу, + 9) 
ит И | 


Исходя из доказанного ранее получаем: 
__" печи, +0)  зш(м(+у,), если т= Ти/Е 
(= Ш —_—_—_еЕ Г 
О и+ 1 0, если у Ти/т; 
7. 0-4+ АИ. 
21 
То есть функция 1; (Г) почти наверное будет равна у (1), при 


Т = ТИ, в остальных случаях некоторые гармонические составляющие 
обратятся в ноль. Но тогда длят Ти 


р,@)- Е Фо+соуш: 


1 9 Е. Ге ь 
= Рош: = Роба - сай © 
Г. Г. Ро 
где С, (®) = але. г эт чу в р„(т) - значение средней за пе- 


риод мощности при Т = Ти. 


Я Т 
Учитывая, что Е а (КС. (Е]4Ё= 0, тогда 


0 
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Т Т 
9 м : си -р@)+ т с. 


0 0 0 


Е" 


Т 


2 
Так как — С. (Е) @Ё> 0, теорема доказана. 

То 

Учитывая, что период интегрирования Т’нам неизвестен, то для 
получения средней за период мощности можно использовать Т’[] Т’, так 
как 


т“ Ь 
вит. 7. ш=р(Т). 


Пример использования алгоритма для диагностирования погреш- 
ности поперечного сечения при токарной обработке. Предваритель- 
ные исследования спектрального состава сигнала виброакустической эмис- 
сии показали, что при токарной обработке стали 20Х инструментом Т15Кб с 
частотой вращения шпинделя 0 = 10 Гц спектр сигнала в частотном 
диапазоне 100-10 кГц имеет две явно выраженные частоты: @, = 2,5 кГци 


0. = 800 Гц. Рассмотрим возможность восстановления периодизированно- 


го с частотой 0 сигнала в наблюдаемом сигнале виброакустической эмис- 
сии (рис.1). 
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Рис.1. Пример реализации сигнала виброакустической эмиссии 


На рис.2,а четко выражены всплески при периоде усреднения 
Т. = 1.25 10°си Ти, где и\ .А также при Т, = 0.4 10°си Ти, 


9 ь 
где и\ .На рис.2,6 при Т; = 10`°с наблюдается максимальный всплеск, 


обусловленный равенством Г.И. = ТИ,, где и.,пМ . На рис.3,4 и 5 
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представлены результаты усреднения соответственно с периодами 
Т= 1.25106; =0410сиТ=10 
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Рис.2. Зависимость средней мощности восстановленного сигнала 
от периода усреднения для функции скорости движения суппорта 
вдоль станины станка 
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Рис.3. Результат усреднения с периодом 1 = 1.25 10`°с 
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Рис.4. Результат усреднения с периодом Т, = 0.4 10`с 


1.6 


х 19 





2.5 




















































































































1 1 | 
0.005 0.01 0.015 


Рис.5. Результат усреднения с периодом тя = 
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Приведенные результаты показывают, что с помощью предложен- 
ного алгоритма удается эффективно определять информацию о периодиче- 
ски повторяющихся с периодом вращения детали возмущениях. 

Выводы: 1. Предложенный алгоритм позволяет определять периодически 
повторяющиеся изменения в сигнале виброакустической эмиссии, который 
имеет следующие прикладные значения: 

- для целей обеспечения качества обработки деталей на металло- 
режущих станках; 

- для целей виброакустической диагностики развивающихся по про- 
странству дефектов в роторных системах (например, турбинах). 

2. Аналитически доказана и экспериментально подтверждена схо- 
димость представленного алгоритма. 
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